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The properties of surface plasmons localized at the interface between graphene and kerr-type nonlinear 
medium in three dimensions are investigated. Compared with surface plasmons at the surface of metal, 
with the inevitable nonlinear refractive effect, the confinement of plasmon can be improved to three times 
than graphene plasmons without nonlinear contribution, but also with almost the same relative propagation 
length. Moreover, the dispersion relation and propation distance of graphene plasmons can be easily con- 
trolled by changing the fermi energy, temperature and relaxation time of graphene. Our results suggest a 
simple but useful potential application for precise nonlinear material sensor using graphene plasmons. 



I. INTRODUCTION 

Graphene with unique optical and electronic 
properties''^ becomes ideal for a number of pho- 
tonic application.'^"''' In particular, graphene plasmons 
(GPs),^'^ which is the collective charge density wave cou- 
pled to electromagnetic field localized between graphene, 
its cover and substrate, have become one of the research 
hotspots.® Compared with surface plasmons in metal, 
CPs provide a conveniently tunable platform for strong 
light-matter interactions, much due to their enormous 
confinement and relatively long propagation length. ^'''^ 
Meanwhile, the induced nonlinear contribution from 
the substrate is always ignored for metallic plasmons 
due to that the nonlinear coefficients are always very 
small. However, this is different for GPs. The strong 
electromagnetic field of GPs (the enhancement factor is 
as large as 10^-10^) makes the nonlinear contribution 
should be included when analyzing the properties of 
graphene plasmons. 

As a truly two-dimensional (2D) electronic sys- 
tem, graphene is a monolayer of carbon atoms closed 
packed in a honeycomb crystal. Previous work verified 
that graphene system can bind surface plasmons with 
transverse-magnetic (TM) polarization when Im{cr} > 0, 
where a is the conductivity of graphene.^ From the fres- 
nel coefficient of the graphene system, the dispersion re- 
lation of graphene system is given as/'^ 

. 4-7r cr e 1 
- r-^ + r^— (1) 

where c is the speed of light in vacuum, e is the dielectric 
function of substrate and 77 is the plasmon wave vector 
normalized to the free light momentum w/c. Because 
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of the momentum 77 is usually much larger than unit and 
dielectric function of substrate, so we use the electrostatic 
limit of this expression and the dispersion of the GPs is 
obtained 



However, when considering the kerr-type nonlinear ef- 
fect of the substrate, Eq.(2) is no longer correct with sim- 
ply changing the dielectric function e with the intensity- 
dependent dielectric expression 

erf = eo+a|E|', (3) 

where a is kerr coefficient of the substrate, and eq are 
the linear and total dielectric function of the substrate, 
E is the electric field, respectively. To solve the problem, 
in this paper we reconsider the dispersion relation of GPs 
under the simple system of an infinite planar graphene 
sheet lying on a dielectric medium, but include the Kerr 
effect of the substrate, and the analytic express of the 
dispersion relation of GPs is obtained. Moreover, the 
dependence of the dispersion relation and propagation 
distance of GPs on the chemical potential, temperature 
and relaxation time of graphehe is also clarified. 



II. SYSTEM DESCRIPTION AND THEORY METHOD 

We consider an infinite-large graphene sheet lying on 
a substrate [Fig. 1(a)], the dispersion property of the 
GPs at the surface between graphene sheet and substrate 
can be calculated as follows. Firstly, graphene is phys- 
ically represented by its in-plane complex conductivity 
cr, which can be calculated based on the local random- 
phase-approximation. The contributions from intraband 
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and interband can be separately expressed as 
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where 17 = hu/^, t = T/j.i, with fi and r are the chem- 
ical potential and a finite relaxation time of graphene, 
respectively. Moveover, T is expressed in units of energy. 
And the total conductivity a can be expressed as 



Cr(a;) = CTlntra + CTinto 



(4) 



Next, for simply the nonlinear response of the sub- 
strate is believed to be isotropic, then the electric field 
in this system can be taken the form as E(r,i) = 
E(2;) exp(i/3x — icut), with /3 is the momentum of GPs. 
Substitute this form to Maxwell equations, we can ob- 
tain: 



-rjHy = edEz 



(5) 



also the electromagnetic field in the air layer z < can 
be written as 



Hy = Be< 



(11) 



where k — describes the confinement of the sys- 

tem. From the boundary continuity conditions of and 
Hy at ^ = 0, we can know that: 



A'(0) ^ kB_ 
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(12) 



then substitute this equation into Eq.(7) at f = we can 
learn that: 



Erf 



+(^)^ (13) 

combining Eq.(9), Eq.(12) and Eq.(13), we can get the 
dispersion equation of GPs finally: 



where ^ = fcz,?7 = /?/fc are normalized length and mo- 
mentum to the vacuum wave-vector fc, respectively. The 
prime is differential to ^. By eliminating the electric field 
terms in this equation array we can easily get the mag- 
netic equation: 
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(6) 



then substitute the electric fields from Eq.(5) to Eq.(3), 
one obtains that 



(7) 



next differentiate Eq.(7) and substitute into Eq.(6), then 
multiply td in both sides, one can get the expression: 
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a OL 
integrate over ^ on both sides, we reach the expression: 
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if we set a and e = 1 ,which means £^(0) eo, the 
Eq.(14) will go back to the Eq.(l) as expected. 

For the case of the normalized momentum 77 is much 
larger than unit and the dielectric functions of substrate, 
the wave-vector of GPs can be finally simplified to the 
following expression: 
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where eci(O) = eo + a |Eo| , and o'(cj) takes the form of 
Eq.(4). 



III. NUMERICAL RESULTS AND DISCUSSION 
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with the equations, the magnetic field can be expressed 
by dielectric function and normalized momentum. 

Assuming that the graphene lies at z = 0, and the 
nonlinear medium occupies the region z > Q. In the 
nonlinear substrate region the electromagnetic field can 
be written as 



E 



Hy=A 
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To understand the nonlinear contribution to GPs more 
clearly, we firstly take the substrate as SiC as an exam- 
ple [Fig. 1(a)]. The parameters of graphene and SiC are 
chosen as follows: the chemical potential, temperature 
and electron relaxation time are set as l.OeV, 300K and 
1.0 X lO^^'^s for graphene, respectively. The experimen- 
tal linear dielectric data of SiC is taken from Palik^^, 
while nonlinear coefficient is chosen as 1.51 x lQ~^'^m^ /V'^ 
from Weber^^ and considered to be non-dispersion. Due 
to the field enhancement factor is the order of 10^ in 
normal GPs,^ the peak of electric field can reach E^ = 
1.6 X lO^^T^^/m^ at the graphene interface. The disper- 
sion relation and propagation distance of GPs without 
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and with nonlinear contribution from the substrate can 
be calculated by Eq. (2) and Eq. (15) with these pa- 
rameters, respectively. And the results are shown in Fig. 
1 (b) . We can see that the peak momentum of GPs almost 
can reach 3000 times to free light momentum (thin solid 
line) when the photon energy (in unit of fermi energy) 
is 1.6 cV, and the relative propagation length can reach 
10 times to SPP wavelength (thin dashed line) without 
nonlinear contribution of substrate. However, if the non- 
linear contribution is taken account, the momentum of 
GPs reach lO'' times to free light wavelength (thick solid 
line), which means that the nonlinear effect can improve 
the localization of electromagnetic filed to 3 times. The 
relative propagation distance of GPs with nonlinear con- 
tribution is slightly larger than GPs without nonlinear 
contribution (thick dashed line), which means that the 
nonlinear effect has little effect in relative propagation 
length. 




FIG. 1. (a) Sketch of an infinite planar graphene sheet ly- 
ing on a substrate, (b) Dispersion relations (sofid lines) and 
propagation distances (dashed lines) of graphene plasmons 
with (thick lines) and without (thin lines) nonlinear contri- 
bution from the substrate. The photon energy is SI = hw/ ^l. 
The momentum is eigenvector in units oi uj/c and the propa- 
gation length given by Im(l/fcsp) in unit of SPP wavelength. 
These magnitudes are obtained for the case of graphene in 
^ = l.QeV,T = 1.0 X 10""s and T = 3007^. The nonlinear 
coefficient is chosen as 1.51 x lO'^^m^/l/^, and the peak of 
electric field is set to El = 1.6 x IQ^'^V'^/m^. 

Moreover, From the Eq. (3), one can see that the non- 
linear part of the dielectric constant of the substrate de- 
pendent not only on the nonlinear coefficient, but also on 
the incident light intensity. So it is natural to consider 
how the nonlinear dielectric function [aE"^) affects the 
properties of the GPs. The nonlinear dielectric function 
resolved dispersion and propagation distance of GPs are 
shown in Fig. 2(a) and (b), respectively. One can see that 
the nonlinear dielectric largely enhances the normalized 
momentum while affects little on the relative propaga- 
tion length, which means that it is possible for us to 
enhance localization of GPs by improving the incident 
electric field or changing to another material with large 
nonlinear coefficient. 

In further, compared with metal, the properties of 
graphene can be controlled by lots of means, such as 
through doping, applying a gate voltage or changing the 
temperature. All of these methods can be efficiently 
change the conductivity of graphene according to Eq. 
(4). which open the active regulation of properties of 
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FIG. 2. The dependence of dispersion relation (a) and relative 
propagation length (b) of graphene plasmons on the nonlin- 
ear dielectric function («-E^) in the system Fig. 1(a). The 
parameters of graphene are set as ^ = l.OeV, r = 1.0 x lO^^^s 
and T = 300 A'. 

GPs for designing applications. Therefore, three most 
important physical parameters for describing graphene, 
chemical potential /x, temperature T and relaxation time 
T are separately studies to look how they affect the prop- 
erties of GPs with nonlinear contributions in Fig. 3, Fig. 
4 and Fig. 5, respectively. 

Firstly, the fermi energy resolved dispersion relation 
and propagation length of GPs are shown in Fig. 3(a) 
and Fig. 3(b), respectively. The normalized momentum 
of GPs reaches maximum when the energy approaches 
1.6 eV. And the higher the fermi energy is, the larger 
the normalized momentum and relative propagation are. 
In addition, the propagation distance of GPs in infrared 
region is much longer than GPs in the visible region, 
though they are not so confinement. These results show 
that higher fermi energy is better for GPs mode, and the 
application region is more likely lying in the infrared or 
terahertz region. 
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FIG. 3. Fermi energy of graphene resolved dispersion relation 
(a) and relative propagation distance (b) of graphene plas- 
mons with nonlinear contribution from the substrate in the 
structure shown in Fig. 1(a). All of the data are calculated 
with r = 1 X 10~^^s, T = 300A', and the nonlinear dielectric 
function is chosen as aEl — 24.1608. 

Secondly, the environmental temperature can also af- 
fect the conductivity of graphene. The collision of elec- 
trons in graphene is directly related to the temperature. 
We calculate the temperature resolved dispersion rela- 
tion [Fig. 4(a)] and propagation length [Fig. 4(b)] of GPs. 
From the contour plots, we know that low temperature 
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benefits to the GPs modes. In lower temperature, GPs 
possess more longer propagation distance, which can ex- 
plained by that the energy loss from coUisons of elec- 
trons in graphene is reduced. More important, it should 
be mentioned that the maximum of momentum reaches 
3.73 X 10^ when the temperature nearly 50 K, this means 
that GPs with not only comparable long propagation dis- 
tance but also extremely field enhancement and localiza- 
tion can realize with low temperature in graphene, which 
can find applications in strong fields physics. 
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FIG. 4. The dispersion relation (a) and relative propagation 
distance (b) of graphene plasmons in the structure sketched 
in FIG. 1(a) with different environmental temperatures. The 
parameters of graphene are set as fj, = l.OeV, r = 1 x 10"^''s, 
and the nonlinear dielectric function is aE^ = 24.1608. 

Thirdly, the electron relaxation time r is also impor- 
tant for the property of graphene. It is proportion to the 
fermi energy in ideal condition, but it is not exact correct 
in the real graphene. Then it opens another freedom to 
adjust T to change the property of GPs. As shown in Fig. 
5(b), the relaxation time affects the propagation length 
of GPs which is similar to electrons in metal. The longer 
the relaxation time is, the longer the plasmon propagates. 
Fig. 5(a) tells us that the relaxation time does not con- 
tribute as much to the dispersion as propagation length 
of the graphene SPs. Except for that, the relaxation time 
affects more in the low energy zone than the higher en- 
ergy region. 
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FIG. 5. Relaxation time of graphene resolved dispersion rela- 
tion (a) and relative propagation length (b) of graphene plas- 
mons with the nonlinear contribution from the substrate (SiC) 
in the structure of Fig. 1(a). These magnitudes are obtained 
for the case of graphene in jj, = leV and T — 300K, and the 
nonlinear dielectric function of SiC is set as oEq = 24.1608. 



IV. DISCUSSION AND CONCLUSION 

Usually, the nonlinear contribution of substrate to GPs 
is ignored due to the small nonlinear coefficient of mate- 
rials. However, things have changed for graphene. The 
very large field localization of GPs leads to super-strong 
field enhancement. In this case the nonlinear contri- 
bution of the material becomes an inevitable factor to 
calculate the surface plasmon dispersion relation even if 
the nonlinear coefficient is very small. Alternatively, this 
property makes GPs potential candidate for the applica- 
tion of precise nonlinear material sensor. In addition, the 
control of the fermi energy, temperature, electron relax- 
ation time of graphene also make the application of GPs 
is more easier for sensor. 

To conclude, in this paper, the dispersion relation and 
propagation distance of GPs with nonlinear contribution 
from the substrate are exactly examined. And compared 
with the case without nonlinear contribution, we found 
that the confinement of GPs can be improved to three 
times with almost the same relative propagation distance. 
At the same time, the fermi energy, temperature and the 
relaxation time of graphene also affect the dispersion and 
relative propagation length. One can find that graphene 
with high fermi energy, low temperature and low realtime 
is the best choice for GPs modes, and these modes can 
propagate longer and have more confinement. The super- 
high field enhancement can be used in nanoparticle sensor 
and fluorescence enhancement. 
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